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ABSTRACT 

We present high angular resolution observations toward two massive star 
forming regions IRAS 18264-1152 and IRAS 23151+5912 with the Plateau de 
Bure Interferometer (PdBI) in the SiO (J=2-l) and Hi3co+ (J=l-0) lines and 
at 1.3 mm and 3.4 mm continuum, and with the Very Large Array (VLA) in 
the NH3 (J,K) = (1,1), (2,2) lines. The NH3 (1,1) and (2,2) emission is detected 
toward IRAS 18264-1152 only. For IRAS 18264-1152, the SiO observations reveal 
at least two quasi-perpendicular outflows with high collimation factors, and the 
most dominant feature is a redshifted jet-like outflow with very high velocities 
up to about Af = 60 kms~^ with respect to the systemic velocity. The very- 
high- velocity component (Af = 22 — 60 kms~^) of this outflow is spatially offset 
from its high- velocity {Av = 3 — 21 kms~^) component. The SiO line profiles 
and position-velocity characteristics of these two components suggest that this 
outflow can be driven by an underlying precessing jet. For IRAS 23151+5912, 
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the bipolar but mainly blueshifted SiO outflow coincides with the inner parts 
of the single-dish CO outfiow. In particular the quasi-parabolic shape of the 
blueshifted outfiow coincides with the near-infrared nebulosity and is consistent 
with entrainment of the gas by an underlying wide-angle wind. The analysis of 
the molecular outfiow data of the two luminous sources further support high-mass 
stars forming via a disk-mediated accretion process similar to low-mass stars. 

Subject headings: ISM: individual (IRAS 18264-1152, IRAS 23151+5912) — ISM: 
jets and outflows — stars: formation 

1. Introduction 

The formation process of massive stars is poorly understood. There are two competing 
views: One view is that high-mass stars form via infall and disk mediated accretion processes 
as a scaled- up version of their low-mass counterparts (e.g., Jijina & Adams 1996; Garay & 
Lizano 1999; Norberg & Maeder 2002; McKee & Tan 2002; Yorke & Sonnhalter 2002; McKee 
& Tan 2003; Keto 2003). The other view argues that at the dense center of evolving massive 
clusters, competitive accretion could become a dominant process. In the extreme case, even 
the coalescence of low- and intermediate- mass (proto) stars may occur to form the most 
massive stellar objects (e.g., Bonnell et al. 1998; Stabler et al. 2000; Bonnell et al. 2004; 
Bally & Zinnecker 2005). 

Massive molecular outflows play a very important role in massive star formation. In 
addition to dissipating excess angular momentum of the infalling material as their low-mass 
counterparts (Shu et al. 1987), they may inject enough energy to sustain the turbulence and 
prevent a decrease of the virial parameter (Quillen et al. 2005), whose value may determine 
whether competitive accretion is effective in cluster forming clouds (Krumholz et al. 2005a). 
Outflows also lead to significant anisotropy in the distribution of the circumstellar material, 
which greatly reduces the radiation pressure experienced by the infalling material (Krumholz 
et al. 2005b). The coUimation factor of massive molecular outflows can help to discriminate 
between the disk mediated accretion process and the coalescence process, since the former 
requires massive outflows as collimated as their low-mass counterparts and the latter predicts 
outflows to be far less collimated. 

Surveys with single-dish telescopes find that massive molecular outflows are ubiquitous 
toward high-mass star forming regions, and are far more massive and energetic than their 
low-mass counterparts (e.g.. Shepherd & Churchwell 1996; Zhang et al. 2001; Beuther et al. 
2002c; Zhang et al. 2005). While Shepherd & Churchwell (1996) find that massive outflows 
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appear to be less coUimated than their low-mass counterparts, Beuther et al. (2002c) argues 
that this may be mostly an observational artifact caused by the large distances of the target 
sources and the low spatial resolution of most studies. In recent years, there have been several 
high spatial resolution mm interferometric studies of molecular outflows toward massive star 
forming regions. Shepherd & Kurtz (1999) reveal a massive outflow with a wide opening 
angle toward G192. 16-3.82 and suggest its driven mechanism as a strong wide-angle wind; 
Cesaroni et al. (1999) report a highly coUimated SiO jet toward IRAS 20126-^4104, and 
Shepherd et al. (2000) reveal the same, but a larger scale (2pc) outflow in CO in a different 
orientation; Hunter et al. (1999) show a well coUimated SiO jet toward AFGL 5142; Beuther 
et al. (IRAS 05358+3543, 2002a; IRAS 19410+2336, 2003; IRAS 19217+1651, 2004) resolve 
simple outflows previously identified with single-dish observations into multiple and well 
coUimated outflows when observed with interferometers at high spatial resolution. They 
find that their kinematics are similar to those of low-mass counterparts; Shepherd et al. 
(2003) present multiple, overlapping massive outflows driven by at least four protostars in 
the region of W75 N. Their energetics and near-infrared observations suggest that they are 
not likely to be scaled-up versions of jet-driven outflows from low-mass protostars; Gibb et 
al. (2003) show that a precessing jet-driven flow associated with G35.2-0.7N in single-dish 
observations can be explained as at least two overlapping flows; Su et al. (2004) reveal a clear 
bipolar morphology of molecular outflows toward two luminous sources IRAS 21519+5613 
and IRAS 22506+5944; Kumar et al. (2004) unveil at least four outflows toward an ultra- 
compact (UC) Hll region Onsala 1; Sollins et al. (2004) resolve a highly energetic bipolar SiO 
outflow toward UC Hll region G5. 89-0. 39. All of these interferometric observations reveal 
new morphologies of massive outflows that can not be resolved with single-dish observations 
and give a wealth of information about the physical process in the innermost parts of the 
massive star forming regions. 

As the statistics of the high spatial resolution interferometric observations arc still poor 
and some basic issues, such as coUimation factors and outflow driving mechanism, are being 
debated, we carry out a study of massive molecular outflows with the Plateau de Bure 
Interferometer (PdBI) and the Very Larege Array (VLA). Here we use the shock tracer SiO 
(J=2-l) to image outflows toward two massive star forming regions IRAS 18264-1152 and 
IRAS 23151+5912 (hereafter we refer as 118264 and 123151, respectively). The dense cores 
and ambient gas are mapped at 1.3 and 3.4 mm continuum, and in the NH3(1,1) and (2,2) 
and H^^CO^ (J=l-0) lines. The two sources are part of a larger sample discussed in detail 
by Sridharan et al. (2002) and Beuther et al. (2002b, d). Beuther et al. (2002c) present 
bipolar CO outflows toward the two sources with the IRAM 30m telescope observations. 
Toward 123151, Weigelt et al. (2006) reveal a cone-hke nebulosity, which coincides with 
the blueshifted CO outflow, in the near-infrared K' band using the bispectrum speckle 
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interferometry method. The kinematic distance of 123151 derived from the CS velocity is 
5.7 kpc. The kinematic distance of 118264, suffering from an ambiguity, is either 3.5 or 12.5 
kpc (Sridharan et al. 2002). By associating the region via the near- and mid-infrared surveys 
2MASS and MSX on larger scales with sources of known distance, the ambiguity was solved 
to be at the near distance (S. Bontemps, 2003, priv com.). On the other hand, toward the line 
of sight of 118264, Sewilo et al. (2004) solve the kinematic distance ambiguity of an Hll region 
with velocity of 50.9 kms~^ to be at the far distance of 12.0 kpc. For the molecular cloud 
with velocity of 43.6 kms"-*^, Sewilo et al. consider it is associated with the Hll region of 50.9 
kms""^ and adopt its far distance of 12.4 kpc. Since the distance determinations above both 
use indirect methods and still give ambiguous results, we give physical parameter estimations 
corresponding to both the far and near distances in the following sections. Based on the 
High Resolution (HIRES) IRAS database, Sridharan et al. (2002) estimate the luminosities 
of the two sources to be IO^Lq for 123151 and 10'^Lq and 1.2 x IO^Lq respectively for the 
near and far distances for 118264. As additional evidence of massive star formation and 
outflows. Class II CH3OH and H2O maser emission toward 118264 and H2O maser emission 
toward 123151 is observed (Beuther et al. 2002d). Sridharan et al. (2002) detect no emission 
at 3.6 cm down to 1 mJy for both sources, while Zapata et al. (2006) detect 1.3 cm and 3.6 
cm emission toward 118264 with a better sensitivity. After a description of observations in 
Sect. 2, we show results in Sect. 3. Discussion is presented in Sect. 4, and we conclude in 
Sect. 5. 



2. Observations 

The observations of 118264 and 123151 were carried out during August 2003 to November 
2003 using the PdBI ^. At the 3mm wave band, the spectral correlator was set to sample the 
SiO (v=0, J=2-l) and H^3C0+ (J=l-0) lines with a bandwidth of 40 MHz. In addition, two 
wide bands of 320 MHz were set to observe the continuum emission, and to cover the SiO and 
H^'^CO^ lines as well. Taking advantage of the dual frequency operation, we set the correla- 
tor to the CN (J=2-l) line using a bandwidth of 40 MHz, and two 320 MHz wide bands at the 
1mm band for continuum. The CN (2-1) emission was not detected. Except for observations 
toward 123151 on 2003 December 9, all the tracks were carried out in weather conditions ade- 
quate for the 3mm band only. 118264 was covered with a 7-field mosaic with the pointing cen- 
ters at {ARA, ADec) = (5", 3".5; 28", 3".5; -18", 3"5; 16"5, 23"5; -6"5, 23"5; 16"5, -16"5; 
-6".5, -16".5) with respect to the reference center at RA(J2000) = 18:29:14.3 and Dec(J2000) 



^Bascd on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported 
by INSU/CNRS (France), MPG (Germany) and IGN (Spain). 
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= -11:50:26. 123151 was observed with a single field with the pointing center at RA (J2000) = 
23:17:21.0 and Dec (J2000) = 59:28:49.00. We used MWC349 as the primary flux calibrator, 
and 3C273/3C345 as the bandpass calibrators. The time dependence gain was monitored 
by observing 1741-038, 1830-210, and 2200-1-420. The visibility data were calibrated and 
imaged using the standard procedure in GILDAS package. The systemic velocities are 43.6 
kms~^ for 118264 and -54.4 kms~^ for 123151. The nominal spectral resolution was about 
0.4 kms~^ and averaged to about 1 kms~^ for SiO (2-1). The continuum rms in the 3mm 
band was ~0.7 mJy for 118264 and ~0.2 mJy for 123151. 

The observations of the NH3 (1,1) and (2,2) lines were carried out on 1997 December 13 
for 123151 in the D configuration, and on 2001 July 23 for 118264 in the C configuration of 
the VLA ^. In both observations, we used the correlator mode 4 which provided 3.13 MHz 
bandwidth for both the left and right polarizations of the (1,1) and (2,2) lines, respectively. 
The spectral resolutions of the observations was 48 KHz, or 0.6 kms~^ at the NH3 fine fre- 
quencies. The pointing center of the observations was RA(J2000) = 18:29:14.31, Dec(J2000) 
= -11:50:25.6 for 118264 and RA(J2000) = 23:17:21.10, Dec(J2000) = 59:28:48.6 for 123151. 
The on-source integration per source was about 1 hour. The visibility data were calibrated 
and imaged in the AlPS package. The left and right polarizations were averaged during 
imaging to reduce the noise level in the data. The synthesized beam size was about 4" for 
the 123151 images, and 2".8 x 2".0 for 118264 images. The rms noise level is 7 mJy per 0.6 
kms~^ channel. 



3. Results 

For both sources, we detected strong emission in SiO (J=2-l) and H^^CO^ (J=l-0). 
118264 has detectable ammonia emission at an rms of 7 mJy per 0.6 kms~^ channel, and 
123151 was not detected in NH3 at an rms of 8 mJy per 0.6 kms~^ channel. 

3.1. Millimeter Continuum sources 

Fig. 1 presents the 3.4 mm and 1.3 mm continuum emission of the two sources. For 
118264, the 1.3 mm and 3.4 mm continuum emission resolve two peaks with the stronger 
one in the west. While for 123151, the continuum emission remains singly peaked even at a 



^The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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resolution of 1".3 x 0".85 at 1.3 mm. For both sources, the peak positions at 1.3 mm and 3.4 
mm coincide well with each other. In 123151, a fan-shaped structure opening to the east can 
be roughly identified in the 3.4 mm continuum emission. In the 1.3 mm continuum emission, 
there is a discontinuous arch structure at the lowest contour level which goes from a little bit 
northeast to the southeast of the mm continuum peak. We will further discuss these features 
in Sect. 3.2. Because the weather conditions were only adequate for the 3 mm wave band 
during the observing seasons, we use the 3.4 mm continuum data for quantitative analysis. 
Assuming that the 3.4 mm continuum is mainly produced by optically thin dust emission, we 
can calculate the masses of the dense cores following the relation Mgas+dust — Fi,D'^ / Bj,{Td)Ki, 
(Hildebrand 1983), where F^, is the flux density of the dust emission, D is the distance to 
the source, and B^, is the Plank function at a dust temperature of T^. Sridharan et al. (2002) 
derive dust temperatures of 35 K for 118264 and 68 K for 123151 by graybody fits to the IRAS 
and mm data. Here the dust opacity per gram is taken to be = 0.1(z//1.2THz)^cm^g^^ 
(Hildebrand 1983), where the opacity index f3 is set to be 1.5. The results of the calculations 
{Mcore) are listed in Table. 1. The uncertainty of this estimation mainly comes from the 
determinations of /3 and T^. The masses will decrease by a factor of 4 if /? = 1, and increase 
by a factor of 2 if decreases to Ta/2. The integrated flux at 3.4 mm is 0.13 Jy for 118264 
and 0.029 Jy for 123151. When compared with the 1.2 mm single-dish flux (Beuther et al. 
2002a), the total flux from PdBI amounts to 63% for 118264 and 56% for 123151 of the 
single-dish flux extrapolated from 1.2 mm using S{i') o^v^^^ with = 1.5. Some extended 
emission is not recovered by the interferometer. 



3.2. SiO outflows 

Fig. 2 presents the channel maps in SiO (2-1) in 118264, where the velocity resolution 
is smoothed to 5 kms~^. The SiO emission mostly appears in the redshifted channels in 
Fig. 2, and only the 38.6 kms~^ channel shows prominent blueshifted emission. The most 
remarkable feature in the channel maps is an elongated structure in the southeast. This 
redshifted emission has very high velocities up to Aw ~60 kms~^ with respect to the systemic 
velocity (vlsr) 43.6 kms~^. We have examined the data from the 320 MHz band and found 
that there is no detectable SiO emission beyond vlsr of 110 kms~^. The integrated blue- 
and redshifted SiO emission is shown in Fig. 3a, where the single-dish bipolar CO outflow 
is resolved into two quasi-perpendicular outflows: The southeast to northwest (SE-NW) 
outflow and the northeast (NE) outflow. Both outflows seem to originate from the western 
peak of the mm continuum. Along the SE-NW outflow, both red- and blueshifted emission 
can be found, which is a typical feature for expanding bow shocks near the plane of the 
sky. From its alignment with the SE jet-like outflow, the bipolar emission in the northwest 
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seems to be the NW lobe of the SE-NW outflow. But it is also possible that this feature 
is due to another low-mass (proto)star whose mass is below our detection limit. The two 
outflows are both well collimated with overall coUimation factors ~ 3 for the NE outflow and 
~ 4 for the SE-NW outflow. The estimated collimation factors should be the lower limits 
considering the unknown inclination angles and the spatial resolution limit. In Fig. 3b, the 
redshifted SiO emission is shown in two velocity ranges, i.e. the high-velocity component 
(HC) with velocities Av = 3 — 21 kms~^, and the very-high- velocity component (VHC) 
with velocities Av — 22 — 60 kms"^. The definition of the velocity ranges here is based 
on the morphological changes in the SiO emission in the velocity channel maps and the 
characteristics in the position-velocity and mass-velocity diagrams that will be discussed 
below. The redshifted SE outflow in the VHC is shifted toward north and its remote peak 
is closer to the driving source than that in the HC. We will discuss this jet-like outflow in 
detail in Sect. 4.3. Thanks to the continuous velocity structures in this dominant outflow, 
we derive the dynamical time following 

tdyn - ^ , 
'-^'^max 

where Lfi^w is the length of the jet- like SE outflow and Av^ax is the maximum velocity of the 
SE outflow. Assuming the optically thin thermal SiO (2-1) emission in local thermodynamic 
equilibrium (LTE), we estimate the gas mass in the outflow according to 

Mout = Nsio [^^] /Ug mn^d^ An, 

where fid is the permanent dipolc moment, AQ is the FWHM of the synthesized beam, 
fig = 1.36 is the mean atomic weight of the gas, 'mH2 is the mass of a hydrogen molecule, 
d is the distance of the source, S^, is the flux density, and other symbols have common 
meanings. The mean excitation temperature T^x is estimated to be 24 K from the NH3(1,1), 
(2,2) emission. It should be noted that [^^] has large uncertainty since the SiO abundance 
can be greatly enhanced by shocks as a result of grain destruction leading to Si injection 
into the gas phase (Scab & ShuU 1983). In contrast to the typical SiO abundance 10^^^ 
to 10"^^ in dark clouds (Ziurys et al. 1989), Mikami et al. (1992) and Zhang et al. (1995) 
detect the SiO enhancement of four to five orders of magnitude toward the LI 157 outflow. 
Hirano et al. (2001) estimate the SiO abundance to be 10^^° to 10^*^ toward the multiple 
outflows in IRAS 16293-2422. We adopt an SiO abundance of 10~^, which corresponds to an 
enhancement of three to four orders of magnitude with respect to that in dark clouds. Then 
the outflow mass can be determined. Consequently we get the outflow rate according to 
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The derived results are listed in Table 1. 

The channel maps with 1 kms~^ resolution in SiO in 123151 are presented in Fig. 4. With 
respect to the systemic velocity -54.4 kms~^ (vlsr): the SiO emission is mainly buleshifted 
and extends to the velocity A^; ~ —14 kms~^. The integrated emission is shown in Fig. 5. 
In addition to the emission close to the mm continuum source, there is a weaker redshifted 
feature to the south. This is maybe an independent outflow from a low-mass (proto)star 
whose mass is below our detection limit. The bipolar SiO outflow seems to emanate from 
the mm continuum peak and coincides with the respective inner parts of the CO outflow 
from the single-dish observations (see Fig. 5a). The dominant buleshifted SiO outflow has 
a quasi-parabolic shape with the mm peak at its tip and a cavity in the center. This 
structure can also be identified in the channel maps, especially in the -56.4 kms~^ channel. 
Using the bispectrum speckle interferometry method, Weigelt et al. (2006) reveal a cone-like 
nebulosity in the near-infrared K' band with a bright near-infrared point source at the tip. 
This nebulosity coincides well with the quasi-parabolic SiO outflow here. In Figs. 5b and 5c 
the quasi-parabolic SiO oTitflow also coincides with the fan-shaped structure in the 3.4 mm 
continuum and the arch structure in the 1.3 mm continuum. The arch at 1.3 mm is detected 
only at 2 — 3o" level. However, the agreement with the 3.4 mm and the near-infrared image 
suggests that the structure is real. Similar to 118264, we derive the characteristic parameters 
of the SiO outflow in 123151, where Tex is assumed to be 30 K. The results are given in Table 1. 



3.3. H^3CO+ condensations 

Figs. 6a and 6b present the H^'^CO^ channel maps with 0.4 kms~^ resolution for 118264 
and 123151, respectively. In both sources, the H^^CO^ emission is limited to velocity chan- 
nels around the systemic velocity, which confirms that this line emission mostly traces the 
dense ambient gas. For 118264, the peaks of the H^^CO"*" emission in the slightly blue- 
and redshifted channels approximately coincide with the western and eastern peaks of the 
mm continuum, respectively. Figs. 7a and 7b present the integrated emission overlaid with 
the SiO outfiows for the two sources. For 118264, the double peaks in the mm continuum 
are roughly resolved in the integrated H^^CO+ emission. In addition to the condensations 
around the mm continuum peaks, the H^^CO'^ emission also shows extensions coincident 
with the SiO outflows. In particular the southeast elongated feature in H^'^CO^ coincides 
well with the jet-like SE SiO outflow. For 123151, the major and minor peaks in the inte- 
grated H^'^CO^ emission correlates well with the SiO clumps. The quasi-parabolic shaped 
SiO outflow also has counterparts in H^^CO"*", especially in the -54.0 kms~^ channel and inte- 
grated emission. The correlation between the integrated H^^CO+ emission and SiO emission 
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for the two sources suggests that the H-'^^CO''" gas may be influenced by the outflows. Adopt- 
ing an HCO"*" abundance of 1 x 10~^ (van Dishoeck et al. 1993) and a C to ^^C ratio of 67 
(Langer & Penzias 1990), we can estimate the gas masses of the H^^CO"'" condensations 
with the similar method used in the estimation of the gas masses in the SiO outflows. The 
results (Mdense) are listed in Table 1. Note that the masses derived from the H^^CO^ emis- 
sion are much larger than the core masses derived from the mm continuum. The integrated 
H^^CO''" emission is much more extended than the interferometric mm continuum for both 
sources. Based on the single-dish 1.2 mm continumm observations, Beuther et al. (2002b) 
derive the masses ~21OOM0 and ~28OOOM0 respectively for the near and far distances for 
118264 and ~62OM0 for 123151, which are consistent with the the masses of the H^^CO+ 
condensations here. Thus the interferometric mm dust continuum probably traces the dens- 
est cores, flltering out most emission and hence yielding less masses, whereas the H^^CO"'" 
traces lower density gas and may even be slightly optically thick. In addition, the H^^CO"*" 
emission is obviously affected by the outflows and HCO^ and H^'^CO^ can be enhanced 
by shocks. Toward the bipolar molecular outflow of L1157, Bachiller & Gutierrez (1997) 
detect the abundance of HCO+ to be enhanced by a factor of 26 — 30. J0rgensen et al. 
(2004) reveal an HCO+ abundance of 2.9 xlO"'^ in the NGC 1333 outflow region. Girart et 
al. (2005) derive an averaged HCO'*' abundance of 3.2x10"^ over the core located ahead of 
HH2. Thus the abundance of HCO"^ can be enhanced by a factor of ~ 3 — 30. Assuming the 
H^^CO^ to HCO+ ratio to be constant, the abundance of H^^CO"'" may be underestimated, 
and consequently the mass is overestimated by a factor of ~ 3 — 30. 



3.4. Ammonia emission 

Sridharan et al. (2002) detected the NH3(J,K)=(1,1), (2,2) inversion lines toward 118264 
with the Effelsberg 100m telescope. The NH3 emission toward 123151 was also detected 
with the Effelsberg telescope, but is about an order of magnitude weaker than for 118264 
(Beuther, priv com.). With the VLA we observed the NH3(1,1) and (2,2) lines toward the 
two regions and detected the NH3(1,1) and (2,2) emission only in 118264. We derive the 
optical depth r(l,l,m) for the main component of the NH3(1,1) line and the rotational 
temperature Trot(2, 2:1,1) (Ho & Townes 1983) at three positions as denoted in Fig. 8a. The 
kinetic temperature Tjti„ can be estimated according to 

"^"^^ ' ■ 1 + (W41.7)M1 + C(2,2^2,l)/C(2,2^1,l)]' 

where C(2,2 2, 1) and C(2,2 1, 1) are the rates of collisional transition between the 
levels (2,2), (2,1) and (1,1) (Walmsley & Ungerechts 1983; Danby et al. 1988). The excitation 
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temperature Tex is determined following 

Ts(l, 1, m) = (Te, - 2.7)(1 - e--(^'i'-)), 

where Tb(1, 1, m) is the brightness temperature of the main component of the NH3(1,1) line. 
Then the gas density n{H2) can be calculated according to 
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1 + 



kin 



where A and C are the Einstein A coefficient and the collision rate, respectively (Ho & 
Townes 1983). The results derived from the ammonia emission are listed in Table 2. The 
derived mean rotation temperature 26 K is a little higher than that of 18 K derived by 
Sridharan et al. (2002), and the mean kinetic temperature 36 K is approximately equal to 
the dust temperature (Sridharan et al. 2002). In Fig. 8, there are peaks in NH3(1,1) and 
(2,2) approximately coincident the mm continuum peaks. It seems that NH3(1,1) and (2,2) 
lines trace a more extended envelope than the mm continuum. The derived gas density is 
~ 1.4 ± 0.6 X 10^ cm"^ and the intrinsic hne width (FWHM) is ~2.5 kms~^ 



4. Discussion 

4.1. Outflow Energetics 

In 118264, the SiO emission reveals two quasi-perpendicular molecular outflows emanat- 
ing from the western mm peak. The mass estimated from the mm continuum is as high as 
570 Mq and 7300 Mq respectively for the near and far distances even with a fraction of the 
flux being filtered out by the interferometer. Since the western peak is much stronger than 
the eastern one, most of the mass is attributed to the western peak. To check the missing 
short spacings, the convolved PdBI SiO data have been compared with the single-dish SiO 
data observed with the IRAM 30m telescope (Beuther, priv com.). We find that at channels 
around the systemic velocity 43.6 kms~^, less than 34% of the SiO flux is filtered out by the 
interferometer. While at channels above vlsr of 50 kms~^, no missing short spacing problem 
occurs. The outflow mass for 118264 in Table 1 is derived from channels of 33.6 — 40.6 kms~^ 
(Mhiue) and 46.6 — 103.6 kms"^ (Mred), which are both 3 kms~^ apart from the systemic 
velocity. So the derived outflow mass for 118264 does not suffer much from missing short 
spacings, and the results are comparable to that derived from the single-dish CO observa- 
tions (Beuther et al. 2002c). Prom the SiO channel maps and the integrated emission Mred is 
mostly due to the jet-hke SE outflow, and gives mass outflow rates of 3.4xlO~^M0 yr'^ and 
1.2x10~^Mq yr~^ for the near and far distances, respectively. With the assumption of the 
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momentum conservation between the observed outflow and the internal jet, and adopting a 
typical jet velocity of 500 kms~^, the mass loss rates caused by the underlying jet or wind 
are respectively about 4.1 x 10~'^MQyr~^ and 1.4 x lO^'^M© ?/r~^. Assuming further a ratio 
between the mass loss rate and accretion rate of approximate 1/3 (Shu et al. 1987; Tomisaka 
1998), we get accretion rates of 1.2 x lO^^M^yr^^ and 4.2 x respectively for 

the near and far distances. Such accretion rates are high enough to overcome the radiation 
pressure of the central (proto)star and form most massive stars (Wolfire & Cassinelli 1987; 
Jijina & Adams 1996; Yorke & Sonnhalter 2002). For 123151, the estimation from the mm 
continuum gives a core of 170 Mq. With the same assumptions as for 118264, the outflow 
rate 2.3 x lO~^M0yr~^ will lead to an accretion rate of 1.9 x 10~^MQyr~^. According 
to the comparison between the convolved PdBI SiO data and that from the IRAM 30m 
single-dish observations (Beuther, priv com.), about 32% — 43% of the fiux is filtered out 
by the interferometer for most velocity channels. Considering the missing flux due to the 
missing short spacings and the outflow mostly being blueshifted, the accretion rate can be 
underestimated. Then the accretion rate for 123151 can also be high enough to overcome the 
radiation pressure and form massive stars (Wolfire & Cassinelli 1987; Jijina & Adams 1996; 
Yorke & Sonnhalter 2002). 



4.2. Mass-velocity diagrams 

Outflows associated with low-mass young stellar objects (YSOs) usually exhibit a "mass 
spectrum" mco{v) oc v~'^ (e.g. Chandler et al. 1996; Lada & Fich 1996). The power law index 
7 is typically ~ 1.7 — 1.8 (Lada & Fich 1996), although the slope often steepens at velocities 
greater than 10 kms~^ from the systemic velocity. From a compilation of 22 sources with 
luminosities ranging from 0.58 Lq to 3 x 10^ Lq, Richer et al. (2000) find that at velocities 
below 10 kms~^ 7 is similar, while at velocities above 10 kms~^, 7 is 3 — 4 for low-mass YSOs 
and 3—8 for luminous YSOs. Toward a sample of 11 objects with luminosities Lboi > lO^L© at 
a distance of 2 kpc. Ridge & Moore (2001) do not find a correlation between mass-spectrum 
slope and bolometric luminosity nor a clear separation between 7s measured for the high 
velocity and low velocity emission. Su et al. (2004) also study the mass-velocity diagrams 
of outflows of two luminous YSOs and find a change in slope at a break-point velocity of 10 
kms"-*^. They suggest that the high velocity {\A.v\ > 10 kms~^) gas may drive the low velocity 
(|Ai>| < 10 kms"-*^) gas. With the assumption of the optically thin SiO thermal emission 
in LTE, we derive the mass-velocity diagrams for the SiO outflows for the two sources (see 
Fig. 9). Since the SiO emission is dominantly redshifted in 118264 and blueshifted in 123151, 
we derive the mass-velocity relations in the corresponding wings for the two sources. The 
redshifted emission in 118264 is mostly attributed to the SE outflow. So Fig. 9a represents 
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the mass-velocity relationship of the HC of the SE outflow and can be fitted by a broken 
power law with the index steepening from 71 = 0.36 ±0.05 to 72 = 1.5 ±0.3 at Av = 10 
kms~^. The mass distribution in the velocity channels of the VHC of the SE outflow has 
large scatters and can hardly be described by a linear fit. Similarly, the blueshifted lobe 
of the SiO outflow in 123151 is also fitted by a broken power law with 71 = 0.9 ±0.1 for 
lA^I < 10 kms~^ and 72 = 3±2 for |Af |>10 kms~^. From three dimensional simulations of 
a dense molecular jet penetrating a dense molecular medium, Smith et al. (1997) predict the 
change in slope at high velocities due to a jet-bow shear layer consisting of molecules which 
survived the jet terminal shock. Models of material being accelerated by jet-driven bow 
shocks (Downes & Ray 1999) have reproduced a power law relationship between mass and 
velocity with 7 increasing with decreasing molecular abundance in the jet. However, Downes 
& Ray (1999) give an upper limit of 3.75 for 7, which is much less than the observed values. 
Therefore the physical origin is still unclear for the broken power law of the mass-velocity 
diagrams. 



4.3. Kinematics 



The SiO in the gas phase of molecular outflows can be produced through the sputtering 
of Si-bearing material in grains, where the sputtering is driven by neutral particle impact on 
charged grains in shocks (Schike et al. 1997). In Fig. 7a, the jet-like SE outflow in 118264 
coincides well with the H^'^CO"'" structure extending to the southeast, suggesting the SiO 
emission in this region may arise from the interaction between the shocks driven by the jet 
or wind from the central (proto)star and the dense ambient gas clump. Such an interaction 
mainly exists in the redshifted lobe since the SiO emission is dominantly redshifted. In 
Fig. 10a, the SiO line at the downstream peak of the HC of the SE outflow (denoted as 
the southern cross in Fig. 3b) shows a profile with a steep decrease toward the systemic 
velocity and a gradual redshifted wing, which suggests the SiO enhancement arising from 
the quiescent material accelerated by the shock when the jet or wind impinges into the dense 
ambient gas. We have checked the corresponding single-dish SiO spectrum and find that such 
a profile also exists. So the steep decrease toward the systemic velocity can not be due to 
the missing short spacings. This kind of SiO profiles and interactions have been observed 
toward several low-mass outflows (e.g. Zhang et al. 1995; Hirano et al. 2001). The fine proflle 
at the downstream peak of the VHC (denoted as the northern cross in Fig. 3b) shows a 
very broad plateau in the redshifted wing with velocities up to about 60 kms~^ from the 
systemic velocity (sec Fig. 10b). Such a broad plateau is consistent with the flattened mass 
spectrum of the VHC in the mass-velocity diagram. For 123151, there is a good correlation 
between the SiO clumps and the integrated H^^CO+ peaks. The SiO line at the peak of the 
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blueshifted lobe (denoted as a cross in Fig. 5) also shows a profile of a steep drop toward 
the systemic velocity and a gradual decline in the blueshifted wing (sec Fig. 10c), and this 
profile is confirmed in the single-dish data. So the SiO emission in 123151 can also arise 
from the interaction between the jet or wind driven by the central (proto)star and the dense 
ambient gas. 

Toward several low luminosity sources (e.g., L1448, Guilloteau et al. (1992); NGC 
2264G, Lada & Fich (1996)), it is observed that the flow velocity increases with distance from 
the driving source. On the contrary, Zhang et al. (2000) find the SiO outfiow velocities in 
LI 157 decreasing toward the clumps farther away from the driving source. Also considering 
the different orientations of the different pairs of SiO clumps, they explain this in terms of 
multiple shocks driven by a processing and episodic jet. The downstream peak of the HC of 
the SE outflow in 118264 is farther away from the driving source than that of the VHC, which 
is similar to the SiO outflow in LI 157. And the VHC of this outflow is northward offset from 
the HC. We suggest that the SE outflow is driven by a precessing jet. In this scenario, the 
lower velocities of the HC as compared with the VHC imply a deceleration of the underlying 
jet after travelling a larger distance; The VHC is younger and driven by newly formed shocks 
when the underlying jet is precessing to the north, and its flattened mass spectrum and the 
very broad line width suggest it does not experience signiflcant deceleration. In addition, 
the VHC of the northwest SiO clump is also slightly shifted to the south of the HC, and 
both the VHC and HC align in line with the VHC and HC of the SE outflow, respectively. 
This also suggests the precessing scenario. Most recently, the VLA observations by Zapata 
et al. (2006) resolve the western mm peak in 118264 into triple sources at 1.3 cm and 7 mm. 
One of the three sources shows a slightly rising spectral index between 3.6 cm and 1.3 cm 
and a flux density at 7 mm larger than the value extrapolated from the 3.6 cm and 1.3 cm 
observations. Zapata et al. (2006) suggest this source to be a combination of a thermal jet 
and disk. This is consistent with our interpretation that the SE outflow being driven by a 
precessing jet. 

In case of 123151, the SiO outflow in 123151 traces the inner parts of the large bipolar 
CO outflow despite the additional feature to the south. As described above, the blueshifted 
outflow is revealed as a quasi-parabolic structure, which is coincident with the similar struc- 
tures in the 3.4 mm continuum and H^^CO"^ line emission and also coincident with the near- 
infrared cone-hke nebulosity by Weigelt et al. (2006). Apparently, this can be interpreted 
as the molecular outflow being entrained by a wide-angle wind, where the low-intensity SiO 
emission in the central region of the cone is due to the cavity cleared by the wind while the 
enhanced SiO emission at the base of the outflow is caused by the interaction between the 
wind and the ambient material in the wall of the cavity. Weigelt et al. (2006) also reveal a 
bright point source at the tip of the cone-like nebulosity. As the point source locating at the 
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center of the single-dish 1.2 mm peak, Weigelt et al. (2006) suggest it to be the protostar 
driving an outflow traced by the nebulosity. The astrometry of the point source is (RA, 
DEC)(J2000) = (23:17:21.02, 59:28:48.0) from the comparison of the speckle image with the 
2MASS images and 2MASS All-Sky Catalog of Point Sources. The position of the mm peak 
of our PdBl observations is (RA, DEC)(J2000) = (23:17:20.88, 59:28:47.7). Considering the 
2MASS positional accuracy of 1", the PdBI mm peak can be the same as the near-infrared 
point source. But it is also possible that the offset is real and the very bright near-infrared 
point source is too evolved to be a protostar and the driving source of the outflow corresponds 
to the interferometric mm peak which is not resolved in the single-dish 1.2 mm continuum 
with an 11" resolution. 

The position- velocity (PV) diagrams of the SiO outflows along the axes marked in Fig. 
3a and Fig. 5a are shown in Fig. 11. For 118264, in addition to the feature around the 
driving source with velocities near the systemic velocity, there are mainly two components 
in the PV diagram: the component with velocities increasing linearly with the distance from 
the driving source; the component farther away from the driving source with maximum 
velocities decreasing with the increasing distance. Compared with the channel maps in Fig. 
2, we can find that these two components correspond to the VHC and HC of the SE outfiow 
in Fig. 3b. There is a pair of faint feature centered at about —27", which corresponds to the 
NW lobe of the SE-NW outflow. The velocity structure of the peak emission of the VHC 
follows the so called "Hubble law" and can be explained in terms of bow-shock acceleration. 
This is compatible with the interpretation that the VHC is excited by the newly formed 
shocks without undergoing signiflcant deceleration. If the broken power law is caused when 
the low velocity (|At;| < 10 kms~^) gas comprises ambient gas entrained by the high velocity 
(|Av| > 10 kms^^) gas (Su et al. 2004), the high velocity gas should inject momentum into 
the low velocity gas and decelerate when it travels downstream. This scenario is in agreement 
with the velocity structure of the HC in the PV diagram of 118264, where the maximum 
velocities decrease with the distance from the driving source. For 123151, the PV diagram 
roughly follows the Hubble law. At the position about 7" from the driving source, emission 
features have a relatively wide range of velocities. 



5. Conclusions 

We carry out a study of SiO outflows toward the high-mass star forming regions 118264 
and 123151. According to the H^'^CO^ observations and SiO line profiles, the SiO emission 
in these two sources can be caused by the interaction between the jet or wind from the 
central young stars and the dense ambient gas. The mass-velocity relation of the HC of the 
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SE outflow in 118264 and the blueshifted outflows in 123151 also can be fltted by a broken 
power law with the slopes steeping at about 10 kms~^ with respective to the systemic velocity 
With the PV characteristic of the terminal velocities decreasing with the distance from the 
driving source, the low velocity {\Av\ < 10 kms^^) gas probably comprises the ambient gas 
entrained by the high velocity (lA^I > 10 kms~^) gas. 

Toward 118264, at least two quasi-perpendicular outflows with high coUimation factors 
(~ 3 — 4) are resolved. The VHC of the SE outflow has very high velocities up to A?; 60 
kms~^. It is north offset from the HC, and its downstream peak is farther away from the 
driving source. Comparing with the outflows in the well studied low-mass source L1157, we 
suggest the SE outflow can be entrained by an underlying precessing jet. The characteristics 
of the HC and VHC of the SE outflow in the mass- velocity and PV diagrams also support 
the processing jet scenario. For 123151, the blueshifted SiO emission traces a quasi-parabolic 
shaped outflow which can be identifled in the 3.4 mm continuum and H^^CO"*" emission and 
coincide with the near-infrared nebulosity in the literature. This outflow can be interpreted 
as the molecular gas entrained by the underlying wide-angle wind. 

The core masses estimated from the 3.4 mm continuum are 570 Mq and 7300 Mq 
respectively for the near and far distances for 118264 and 170 Mq for 123151. With the 

assumption of the momentum conservation between the outflow and the driving agent, the 
estimated outflow rates lead to accretion rates of 1.2 x IO^'^'MqUT^^ and 4.2 x lO^'^M^yr^^ 
respectively for the near and far distances for 118264 and 1.9 x 10^^ MQyr^^ for 123151. 
Taking into account the missing flux due to the the missing short spacings for 123151, the 
accretion rate in this region should be higher. Thus for both sources the accretion rates are 
high enough to overcome the radiation pressure from the central objects and form massive 
stars. 

To summarize, our molecular outflow data derived from high resolution observations 
toward two luminous sources show morphologies and kinematics similar to those of the low- 
mass sources. The outflows can be interpreted by jet or wide-angle wind entrainment models. 
The presented data and analysis further support that massive stars up to IO^Lq (~ 30Mq) 
form via disk mediated accretion processes as low-mass stars. 

We are grateful to F. Geuth and the staff at PdBI for their help in the observations and 
data reduction. K. Q. acknowledges the support of the Grant 10128306 from NSFC. H. B. 
acknowledges flnancial support by the Emmy-Noether-Program of the Deutsche Forschungs- 
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Fig. 1. — Top: the 3.4 mm (left) and 1.3 mm (right) continuum emission in 118264. The 
lowest and spacing contours are 2.1 mJy (3cr) at 3.4 mm and 22.5 mJy (3cr) at 1.3 mm. The 
dashed contours represent the —3a level. Bottom: similar to the top but for 123151. The 
lowest and spacing contours at 3.4 mm are 0.6 mJy (3(7). At 1.3 mm, the contours start at 
4.25 mJy {2.5a) in steps of 5.1 mJy (3a-). The dashed contours represent the —3a level for 
3.4 mm and —2.5a level for 1.3 mm. The stars hereafter mark the 1.3 mm continuum peaks. 
The beams are shown as thick eUipses at the lower left of each panel. 
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Fig. 2. — The SiO channel maps in 118264, with a smoothed velocity resolution of 5 kms~^. 
The contours start at 33 mJy (3(t) in steps of 22 mJy {2a). The dashed contours represent 
the -3(7 level. The numbers in the top left of each panel denote the central velocity of each 
channel. The thick ellipse in the lower left of the first panel shows the beam size. 
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Fig. 3. — The integrated SiO emission in 118264. (a) The blue and red contours show the 
blue- and redshifted SiO emission, respectively. The blueshifted emission is integrated from 
33.6 kms~^ to 40.6 kms~^ with contour levels starting at 0.21 Jy-kms~^ (3cr) in steps of 0.14 
Jy-kms~^ (20"). The redshifted emission is integrated from 46.6 kms^^ to 104.6 kms~^ with 
contour levels starting at 0.6 Jy-kms~^ (3(t) in steps of 0.4 Jy-kms~^ (2cr). The blue and red 
dashed contours represent the -3(T level of the blue- and redshifted emission, respectively. The 
solid and dashed grey contours represent the blue- and redshifted lobes of the CO outflow 
in single-dish observations (Beuther et al. 2002c). The solid straight line denotes the axis 
along which the PV diagram is plotted, (b) The gray and red contours denote the redshifted 
emission integrated from the HC channels and from the VHC channels, respectively. The 
gray contours start at 0.339 Jy-kms~^ (3cr) in steps of 0.226 Jy-kms~^ (2o"), and the red 
contours start at 0.486 Jy-kms~^ (3cr) in steps of 0.324 Jy-kms^^ (2o"). The two crosses mark 
the sites at which the SiO spectra are drawn. 
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Fig. 4. — Same as in Fig. 2, but for 123151 with 1 kms ^ resolution. The contours start at 
16.5 mJy {3a) in steps of 11 mJy {2a). 
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Fig. 5. — The blue and red contours represent the blue- and redshifted SiO emission in 
123151, which are integrated from -68.4 kms~^ to -55.4 kms~^ and from -53.4 kms~^ to -50.4 
kms~^, respectively. The blue contours star at 0.115 Jy-kms~^ {5a) in steps of 0.046 Jy-kms~^ 
(2cr) and the red contours start at 0.08 Jy-kms~^ {5a) in steps of 0.032 Jy-kms^^ (2o"). The 
solid and dashed grey contours in (a) describe the blue- and redshifted CO outflow in single- 
dish observations (Beuther et al. 2002c). The solid straight line and the cross mark the axis 
for PV plotting and the site for the SiO spectrum drawing,respectively. The grayscales in 
(b) and (c) represent the 3.4 mm and 1.3 mm continuum emission, respectively. 
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Fig. 6. — (a) The H^^CO^ (1-0) channel maps for 118264. The lowest and spacing contours 
are 78 mJy (Sc). The dashed contours represent the -3o" level. The thick ellipse in the lower 
left of the first panel shows the beam size, (b) Same as for (a), but for 123151. The contours 
start at 21 mJy (Scr) in steps of 14 mJy {2a). 
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Fig. 7. — (a) The grayscale and contours represent the H^^CO'^ (1-0) emission integrated 
from 41.2 to 45.4 kms~^ for 118264, starting at 192 mJy-kms~^ (4o") in steps of 144 mJy-kms~^ 
{3a). The dashed contours represent the -Aa level. The blue and red contours show the SiO 
outflow same as in Fig. 3a but at 5a, 10a, 15a, 20a, 25a levels, (b) Same as for (a), but for 
123151. The grayscale and contours represent the H^^CO"'" (1-0) emission integrated from 
-56.0 to -52.8 kms~^, starting at 36 mJy-kms~^ {3a) in steps of 36 mJy-kms~^. The blue 
and red contours represent the integrated SiO emission same as in Fig. 5 but at 5a, 9a, 13a, 
17a, 21a, 25a levels. 
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Fig. 8. — (a) The NH3(1,1) emission in 118264 integrated from 41.8 kms~^ to 45.4 kms~^. 
The contours start at 20 mJy-kms~^ in steps of 13.4 mJy-kms~^. The pluses denote the 
positions at which the physical parameters in Table. 2 are calculated, (b) Same as for (a), 
but for NH3(2,2) emission integrated from 41.6 kms~^ to 44.6 kms~^. The levels of contours 
are from 19.5 mJy-kms"^ in steps of 13 mJy-kms~^. The dotted contours represent the -3a 
and -6a levels. The grayscale represent the 3.4 mm continuum. The beam sizes are marked 
at the lower left of each panel. 



-28- 




Fig. 9. — (a) Outflow masses calculated for the near distance for 1 kms~^ channels as a 
function of flow velocities for 118264. Note outflow masses for the far distance will give the 
same broken power law. (b) Same as for (a), but for 123151. The sohd lines show the power 
law flts of m.{v) oc v~'^. 
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Fig. 10. — The SiO line profiles. (a)The SiO line profile of 118264 at the position denoted as 
the southern cross in Fig. 3b; (b)Same as for (a), but at the position denoted as the northern 
cross in Fig. 3b; (c)The SiO hne profile of 123151 at the position denoted as a cross in Fig. 
5a. 
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Fig. 11. — (a) The SiO position- velocity diagram for 118264 along the axis denoted as a 
straight line in Fig. 3a. The greyscale and solid contours start at 2a in steps of 2a. The 
dashed contours mark the -2a level, (b) Same as (a), but for 123151 and along the axis 
denoted as a straight line in Fig. 5a. The contour leves start at 2.ba in steps of 2a and the 
dashed contours represent the —2.5(7 level. 
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Table 1. Distances and derived physical parameters for the two sources 



source 


distance 




^dense^ 






Mout" 








(kpc) 




(Mo) 




(Mq) 


(Mo) 


(lO-^MQ^r-i) 


(10* yr) 


118264 


3.5 


570 


3900 


3.2 


17 


20.2 


34 


0.5 




12.5 


7300 


50000 


41 


220 


261 


120 


1.8 


123151 


5.7 


170 


860 


2.2 


0.5 


2.7 


2.3 


1.2 



^masses of the dust-gas cores derived from the 3.4 mm continuum emission 
''masses of the dense ambient gas derived from the H^^CO+ emission 
"^outflow masses derived from the blueshifted SiO emission 
''outflow masses derived from the redshifted SiO emission 

''total outflow masses derived from the sum of the bule- and redshifted outflow masses 
'mass outflow rates derived from M^^^i for 118264 and from Mout for 123151 
^dynamical time scale 
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Table 2. Physical parameters derived from the ammonia emission for 118264. 



position'' 


r(l,l,m) 


T^ot(2,2 : 1,1) 






n{H2) 








(K) 


(K) 


(K) 


(10^ cm-3) 


(kms-l) 


left 


4.0 


24 


23 


30 


1.6 


3.1 


middle 


2.7 


30 


25 


45 


1.0 


1.9 


right 


2.7 


24 


24 


32 


1.6 


2.5 


mean'^ 


3.1 


26 


24 


36 


1.4 


2.5 



"marked as pluses in the left panel of Fig. 8 
'^averaged from the three positions 



